Suzaku Observations of the Radio Galaxy 3C33 
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Abstract 

We present results from a new 100-ks Suzaku observation of the nearby radio galaxy 3C33, and in- 
vestigate the nature of absorption, reflection, and jet production in this source. We model the 2-70 keV 
nuclear continuum with a power law that is absorbed either through one or more layers of pc-scale neutral 
material, or through a modestly ionized pc-scale obscurer. The expected signatures of reflection from a 
neutral accretion disk are absent in 3C 33: there is no evidence of a relativistically blurred Fe Ka emission 
line, and no Compton reflection hump above 10 keV. We discuss the implications of this for the nature of 
jet production in 3C33. 
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1. Overview 

The origin of jets in active galactic nuclei (AGN) is one 
of the most important unsolved problems in extragalac- 
tic astrophysics. While 90% of all AGN (Seyfert galaxies 
and radio-quiet quasars) show little or no jet emission, 
the remaining 10% (the radio-loud AGN and radio-loud 
quasars) launch powerful, relativistic twin jets of parti- 
cles from their cores. Since jets transport a significant 
fraction of the mass-energy liberated during the accre- 
tion process, sometimes out to ~Mpc distances, under- 
standing how they are produced is key to a complete 
picture of accretion and feedback in AGN. 

X-ray observations of the nuclei of radio-loud and 
radio-quiet AGN are essential for establishing the con- 
nection between the accretion flow and jet. Continuum 
X-ray observations of radio-loud AGN have mostly been 
restricted to bright broad-line radio galaxies (BLRGs) 
and quasars, which are oriented relatively close to the 
line of sight with respect to the observer. In these 
sources, unabsorbed non-thermal emission from the jet 
could potentially contaminate the unabsorbed accretion- 
related X-ray spectrum and thus dilute the apparent 
strength of the Compton reflection component. Narrow- 
line radio galaxies (NLRGs), such as 3C 33, on the other 
hand, which are oriented at low to intermediate an- 
gles, have the distinct advantage that (unabsorbed) jet- 
related X-ray emission can be readily spectrally sep- 



arated from (heavily absorbed) accretion-related emis- 
sion, allowing a direct measurement of the strength of 
Compton reflection. However, narrow-line radio galax- 
ies tend to be relatively faint X-ray sources compared to 
broad-line objects. Suzaku observations of these sources, 
therefore, are particularly useful, owing to the high effec- 
tive area of the X-ray Imaging Spectrometer (XIS) and 
Hard X-ray Detector (HXD) instruments. 



Here, we present the results from a 100-ks Suzaku ob- 
servation of one of the brightest NLRGs, 3C 33, the only 
such source known so far to show potential evidence for 
Compton reflection in its 2-10 keV X-ray spectrum. Pre- 
vious observations of 3C 33 with Chandra and XMM- 
Newton (Kraft et al. 2007) showed that its continuum 
spectrum could not be adequately modeled by the combi- 
nation of a heavily obscured power law and a soft power 
law normally fitted to narrow-line FRII sources, due to 
the large residuals present between 2-4 keV. They sug- 
gested that one possibility is that the source is reflection- 
dominated, and that a Compton-reflection continuum is 
present in the 2-10 keV spectrum. However, the limited 
energy range of Chandra and XMM-Newton meant that 
the evidence for reflection in 3C 33 could be unambigu- 
ously claimed. 
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Fig. 1. Suzaku XIS Fl (black), XIS Bl (red), and HXD PIN (green) 
spectra and residuals of 3C33. The model fit is double partially 
covered primary law of photon index T = 1.77, and a Gaussian 
Fe Ka emission line. 

2. Suzaku Observations 

We observed 3C33 with Suzaku on 2007 December 26 
(OBSID 702059010) for a nominal exposure of 100 ks. 
Both the X-ray Imaging Spectrometer (XIS) and Hard 
X-ray Detector (HXD) were operated in their normal 
modes. The source was positioned at the nominal aim- 
point of the HXD instrument. The data were processed 
using v. 2.1.16 of the Suzaku processing pipeline, which 
includes the latest Charge Transfer Inefficiency (CTI) 
correction applied for the XIS. We used the standard 
cleaned events files. We extracted the spectrum of 3C 33 
from the XIS CCDs using a source-centered circle of ra- 
dius 2.5', with background sampled from an adjacent 
region free from any unrelated sources, as well as the 
55 Fe calibration sources at the corners of each detec- 
tor. For the HXD, we added the NXB and CXB spec- 
tra together using the MATHPHA tool. Using the appro- 
priate response file for epoch 4 of Suzaku observations, 
ae_hxd_pinhxnome4_20080129.rsp, we binned the source 
spectrum to 3cr above the total (NXB + CXB) back- 
ground level. 

3. Suzaku Spectral Analysis 

We restricted the energy range for our spectral fits to 2- 
10 keV (XIS) and 15-70 keV (PIN). We fitted to 3C33 
our canonical model for the X-ray spectrum of NLRGs: 
the combination of a heavily absorbed power law (likely 
to be associated with the accretion flow) and a soft, un- 
absorbed power law (Evans et al. 2006). The unabsorbed 
power law is simply a parametrization and we defer the 
discussion of its origin until Section 4.. This model gave 
a relatively poor fit to the spectrum: \ 2 — 150 for 116 
dof, with clear residuals at ~3 keV. The photon index of 
the primary power law T = 1.45 is atypically flat with 
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Fig. 2. Suzaku XIS Fl (black), XIS Bl (red), and HXD PIN (green) 
spectra and residuals of 3C 33 in the energy range 2—70 keV. The 
model consists of a power law modified by a single zone of warm 
absorption, together with a soft, unabsorbed power law, and a 
Gaussian Fe Ka emission line. 

respect to similar radio galaxies studied by Evans et al. 
(2006), and the strong residuals at energies ~3 keV likely 
necessitate a more complex spectral model. We discuss 
such models next. 

3.1. Model I: Complex Neutral Absorption 

We investigated the possibility that an additional layer 
of cold absorption is required to fit the broadband spec- 
trum of 3C 33 adequately. We adopted a double partial- 
covering model, implemented in XSPEC/ISIS as 
phabsxzpcfabs(l)xzpcfabs(2) x (powerlaw+zgauss). We 
included a Gaussian to represent the neutral Fe Ka 
line. This model provides a good fit to the spec- 
trum (% 2 =112 for 114 dof). The best-fitting param- 
eters of this model are a power law of photon index 
T=1.77taH, absorbed by columns N n>1 =(2.1+°- 6 5 ) x 
10 23 cm" 2 (/ 1= 91+»%) and N^H^tl'l) x 10 23 cm" 2 
(/ 2 =80tg%)- Figure 1 shows the Suzaku spectrum and 
residuals for this model. 

We determined the strength of any Compton reflection 
in the spectrum of 3C 33 by replacing the primary power 
law with a PEXRAV component. This neutral reflection 
model provided no improvement in the fit (x 2 = 112 for 
113 dof). The best-fitting value of the reflection fraction 
is R = 0, with a 90%-confidence upper limit of Rj0.5. 

3.2. Model II: Warm Absorption 

An alternative method of modeling the observed spec- 
tral complexity in 3C 33 is to allow the absorber to be 
partially ionized. We modeled the Suzaku spectrum as a 
power law, covered by a single zone of partially photoion- 
ized absorption, plus a second, unabsorbed power law to 
account for the low-energy flux. We included a Gaussian 
to represent the neutral Fe Ka line. This model provided 



an excellent fit to the spectrum (x 2 = 110 for 115 dof), 
with a column density 7Vh~7.7x 10 23 cm~ 2 of modestly 
ionized gas [log(£)~l.l ergs cm s _1 ] as its best-fitting 
absorption parameters. The best-fitting photon index of 
the power law is r = 2. 00 ±0.12. We show the Suzaku 
spectrum, this model, and residuals of x m Figure 2. 
Again, we determined the strength of any Compton re- 
flection in the spectrum of 3C33, and found the best- 
fitting value of the reflection fraction to be R = 0, with 
a 90%-confidence upper limit of RjO.5. 

3.3. Fe Ka Diagnostics 

The strong Fe Ka line detected with Suzaku allows us 
to place constraints on the location and physical state 
of the fluorescing material. The energy of the line core, 
6.385±° 

.023 keV, is consistent with Fe fluorescence from 
neutral or near-neutral (up to ~Fe XVIII) species. The 
unresolved Fe Ka line width and the estimated black 
hole mass of -4 x 10 8 M© (Smith et al. 1990; Bettoni 
et al. 2003) provide a lower limit to the emission radius 
of i; 0.1 pc (^2000 R s ), using Keplerian arguments. Fi- 
nally, the equivalent width of the Fe Ka line is consistent 
with transmission through an absorbing column of rela- 
tive iron abundance Ap c —l.O and column density similar 
to that measured from our spectral fitting (Miyazaki et 
al. 1996). There is no evidence for an additional velocity 
broadened component of the Fe Ka line, which rules out 
the presence of relativistically blurred fluorescent emis- 
sion from the innermost portions of the accretion flow. 
This may imply that a standard, cold Shakura-Sunyaev 
accretion disk is truncated at a large distance from the 
central supcrmassive black hole. We cannot rule out the 
presence of a modestly broadened (R ~ 100i? s ), ionized 
Fe Ka line in 3C 33 but nonetheless it seems clear that 
fluorescent emission from the inner disk is absent. 

4. Soft X-ray Emission: Continuum or Photoionized Gas? 

There has been considerable debate as to whether the 
soft X-ray emission in 3C 33 is dominated by continuum 
emission from the accretion flow or jet (e.g., Kraft et 
al. 2007), or rather originates in a photoionized plasma 
(Torresi et al. 2009). The deep Suzaku spectrum helps 
us to resolve this issue. In Figure 3 we show the co- 
added Suzaku FI and BI spectrum from our observa- 
tion, grouped to a minimum signal-to-noise of 3 and to 
their HWHMs. We also plot the potential transitions 
that would indicate that the the soft-band X-ray emis- 
sion originates in an ionized plasma. For a photoionized 
plasma, these features include H- and He-like species, as 
well as narrow Radiative Recombination Continuua (e.g. 
Evans et al. 2006b). There is tentative evidence for He- 
like transitions of Ne IX that were first reported by Tor- 
resi et al. (2009), and potentially some structure around 
the He-like triplet of Mg XI. Detailed spectral fitting in 



this band, however, is necessary to investigate in detail 
evidence for photoionized gas in the Suzaku spectrum of 
3C33. 

5. Interpretation: Compton Reflection in 3C33 and other 
Radio-Loud AGN 

Our results demonstrate that 3C33 shows no signs of 
Compton reflection from neutral material in the inner 
regions of an accretion disk: there is no reflection hump 
at energies >10 keV, and no evidence for relativistically 
broadened Fe Ka emission. Indeed, the narrow Fe Ka 
line is likely to originate in a much farther region, at 
least 2,000 R s from the black hole. 

These characteristic features of inner-disk reflection, 
which unified models of AGN predict to be prominent, 
are commonly observed in radio-quiet sources (Seyfert 
galaxies and radio-quiet quasars) (Reeves et al. 2006; 
Nandra et al. 2007), but appear to be systematically 
weaker in radio- loud AGN (e.g., Reeves & Turner 2000; 
Evans et al. 2006), as we see in 3C33. 

Several models have been proposed to explain the 
weakness of Compton reflection continuua in radio-loud 
AGN, including the idea that accretion disks are highly 
ionized in high M systems (Ballantyne et al. 2002). In 
this case, the reflection fraction remains high (ii~l), but 
the ionized accretion-disk surface results in a low mea- 
sured R. An alternative interpretation (Garofalo et al. 
2009) is that the black holes in 'high-excitation radio 
galaxies', such as 3C 33, have retrograde spin with re- 
spect to their accretion disks. This pushes the innermost 
stable circular orbit (ISCO) of the accreting material 
outwards, meaning that inner-disk reflection is greatly 
suppressed. Although we cannot distinguish between 
the two models in the case of 3C33, we advertise here 
that we are currently working on a global interpretation 
of the cosmological evolution and implications of black- 
hole spin in radio galaxies (Garofalo, Evans, & Sambruna 
2009). 
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Fig. 3. Co-added Suzaku XIS Fl and Bl spectrum of 3C 33 in the energy range 0.5—1.5 keV. The data have been grouped to a minimum 
signal-to-noise of 3 and to their HWHMs. Overplotted are the lines one would expect to detect in a photoionized plasma, including RRCs, 
and H- and He-like species. 



